Tropical super flies: Integrating Cas9 into Drosophila ananassae and its phenotypic effects.

Ectotherms such as insects are animals whose body temperature largely depends on ambient temperature and temperature variations provide a selection pressure affecting the geographical distribution of these species. However, over the course of evolution, some insect species managed to colonize environments characterized by various temperature ranges. Therefore, insects provide an excellent study system to investigate the basis of adaptation to temperature changes and extremes. We are generally using the vinegar fly Drosophila ananassae as a model system to investigate the genetic basis of cold tolerance. This species has expanded from its tropical ancestral range to more temperate regions resulting in a cosmopolitan, domestic distribution. Previously, we identified candidate genes significantly associated with cold tolerance in this species. We now established molecular genetic tools to assess the function of these genes. Using CRISPR/Cas9 methodology for genome editing and the PiggyBac system, the Cas9 enzyme was successfully integrated into the genome of three fly strains with different levels of cold tolerance. We further report on preliminary findings that the Cas9 integration itself did not have a consistent effect on tolerance to cold. In conclusion, we offer with our study the molecular tools that allow studying stress-related candidate genes in D. ananassae in the future. In addition, we point out and provide guidance on the challenges that come with genome editing in a non-model species.

RelB contributes to the survival, migration and lymphomagenesis of B cells with constitutively active CD40 signaling.

Activation of CD40-signaling contributes to the initiation, progression and drug resistance of B cell lymphomas. We contributed to this knowledge by showing that constitutive CD40-signaling in B cells induces B cell hyperplasia and finally B cell lymphoma development in transgenic mice. CD40 activates, among others, the non-canonical NF-ĸB signaling, which is constitutively activated in several human B cell lymphomas and is therefore presumed to contribute to lymphopathogenesis. This prompted us to study the regulatory role of the non-canonical NF-ĸB transcription factor RelB in lymphomagenesis. To this end, we crossed mice expressing a constitutively active CD40 receptor in B cells with conditional RelB-KO mice. Ablation of RelB attenuated pre-malignant B cell expansion, and resulted in an impaired survival and activation of long-term CD40-stimulated B cells. Furthermore, we found that hyperactivation of non-canonical NF-кB signaling enhances the retention of B cells in the follicles of secondary lymphoid organs. RNA-Seq-analysis revealed that several genes involved in B-cell migration, survival, proliferation and cytokine signaling govern the transcriptional differences modulated by the ablation of RelB in long-term CD40-stimulated B cells. Inactivation of RelB did not abrogate lymphoma development. However, lymphomas occurred with a lower incidence and had a longer latency period. In summary, our data suggest that RelB, although it is not strictly required for malignant transformation, accelerates the lymphomagenesis of long-term CD40-stimulated B cells by regulating genes involved in migration, survival and cytokine signaling.

Natural variation in the transcriptional response of Drosophila melanogaster to oxidative stress.

Broadly distributed species must cope with diverse and changing environmental conditions, including various forms of stress. Cosmopolitan populations of Drosophila melanogaster are more tolerant to oxidative stress than those from the species' ancestral range in sub-Saharan Africa, and the degree of tolerance is associated with an insertion/deletion polymorphism in the 3' untranslated region of the Metallothionein A (MtnA) gene that varies clinally in frequency. We examined oxidative stress tolerance and the transcriptional response to oxidative stress in cosmopolitan and sub-Saharan African populations of D. melanogaster, including paired samples with allelic differences at the MtnA locus. We found that the effect of the MtnA polymorphism on oxidative stress tolerance was dependent on the genomic background, with the deletion allele increasing tolerance only in a northern, temperate population. Genes that were differentially expressed under oxidative stress included MtnA and other metallothioneins, as well as those involved in glutathione metabolism and other genes known to be part of the oxidative stress response or the general stress response. A gene coexpression analysis revealed further genes and pathways that respond to oxidative stress including those involved in additional metabolic processes, autophagy, and apoptosis. There was a significant overlap among the genes induced by oxidative and cold stress, which suggests a shared response pathway to these two stresses. Interestingly, the MtnA deletion was associated with consistent changes in the expression of many genes across all genomic backgrounds, regardless of the expression level of the MtnA gene itself. We hypothesize that this is an indirect effect driven by the loss of microRNA binding sites within the MtnA 3' untranslated region.

Drosophila Evolution over Space and Time (DEST): A New Population Genomics Resource.

Drosophila melanogaster is a leading model in population genetics and genomics, and a growing number of whole-genome data sets from natural populations of this species have been published over the last years. A major challenge is the integration of disparate data sets, often generated using different sequencing technologies and bioinformatic pipelines, which hampers our ability to address questions about the evolution of this species. Here we address these issues by developing a bioinformatics pipeline that maps pooled sequencing (Pool-Seq) reads from D. melanogaster to a hologenome consisting of fly and symbiont genomes and estimates allele frequencies using either a heuristic (PoolSNP) or a probabilistic variant caller (SNAPE-pooled). We use this pipeline to generate the largest data repository of genomic data available for D. melanogaster to date, encompassing 271 previously published and unpublished population samples from over 100 locations in >20 countries on four continents. Several of these locations have been sampled at different seasons across multiple years. This data set, which we call Drosophila Evolution over Space and Time (DEST), is coupled with sampling and environmental metadata. A web-based genome browser and web portal provide easy access to the SNP data set. We further provide guidelines on how to use Pool-Seq data for model-based demographic inference. Our aim is to provide this scalable platform as a community resource which can be easily extended via future efforts for an even more extensive cosmopolitan data set. Our resource will enable population geneticists to analyze spatiotemporal genetic patterns and evolutionary dynamics of D. melanogaster populations in unprecedented detail.

The discovery, distribution, and diversity of DNA viruses associated with Drosophila melanogaster in Europe.

Drosophila melanogaster is an important model for antiviral immunity in arthropods, but very few DNA viruses have been described from the family Drosophilidae. This deficiency limits our opportunity to use natural host-pathogen combinations in experimental studies, and may bias our understanding of the Drosophila virome. Here, we report fourteen DNA viruses detected in a metagenomic analysis of 6668 pool-sequenced Drosophila, sampled from forty-seven European locations between 2014 and 2016. These include three new nudiviruses, a new and divergent entomopoxvirus, a virus related to Leptopilina boulardi filamentous virus, and a virus related to Musca domestica salivary gland hypertrophy virus. We also find an endogenous genomic copy of galbut virus, a double-stranded RNA partitivirus, segregating at very low frequency. Remarkably, we find that Drosophila Vesanto virus, a small DNA virus previously described as a bidnavirus, may be composed of up to twelve segments and thus represent a new lineage of segmented DNA viruses. Two of the DNA viruses, Drosophila Kallithea nudivirus and Drosophila Vesanto virus are relatively common, found in 2 per cent or more of wild flies. The others are rare, with many likely to be represented by a single infected fly. We find that virus prevalence in Europe reflects the prevalence seen in publicly available datasets, with Drosophila Kallithea nudivirus and Drosophila Vesanto virus the only ones commonly detectable in public data from wild-caught flies and large population cages, and the other viruses being rare or absent. These analyses suggest that DNA viruses are at lower prevalence than RNA viruses in D.melanogaster, and may be less likely to persist in laboratory cultures. Our findings go some way to redressing an earlier bias toward RNA virus studies in Drosophila, and lay the foundation needed to harness the power of Drosophila as a model system for the study of DNA viruses.

Genomic Analysis of European Drosophila melanogaster Populations Reveals Longitudinal Structure, Continent-Wide Selection, and Previously Unknown DNA Viruses.

Genetic variation is the fuel of evolution, with standing genetic variation especially important for short-term evolution and local adaptation. To date, studies of spatiotemporal patterns of genetic variation in natural populations have been challenging, as comprehensive sampling is logistically difficult, and sequencing of entire populations costly. Here, we address these issues using a collaborative approach, sequencing 48 pooled population samples from 32 locations, and perform the first continent-wide genomic analysis of genetic variation in European Drosophila melanogaster. Our analyses uncover longitudinal population structure, provide evidence for continent-wide selective sweeps, identify candidate genes for local climate adaptation, and document clines in chromosomal inversion and transposable element frequencies. We also characterize variation among populations in the composition of the fly microbiome, and identify five new DNA viruses in our samples.

Sawfly Genomes Reveal Evolutionary Acquisitions That Fostered the Mega-Radiation of Parasitoid and Eusocial Hymenoptera.

The tremendous diversity of Hymenoptera is commonly attributed to the evolution of parasitoidism in the last common ancestor of parasitoid sawflies (Orussidae) and wasp-waisted Hymenoptera (Apocrita). However, Apocrita and Orussidae differ dramatically in their species richness, indicating that the diversification of Apocrita was promoted by additional traits. These traits have remained elusive due to a paucity of sawfly genome sequences, in particular those of parasitoid sawflies. Here, we present comparative analyses of draft genomes of the primarily phytophagous sawfly Athalia rosae and the parasitoid sawfly Orussus abietinus. Our analyses revealed that the ancestral hymenopteran genome exhibited traits that were previously considered unique to eusocial Apocrita (e.g., low transposable element content and activity) and a wider gene repertoire than previously thought (e.g., genes for CO2 detection). Moreover, we discovered that Apocrita evolved a significantly larger array of odorant receptors than sawflies, which could be relevant to the remarkable diversification of Apocrita by enabling efficient detection and reliable identification of hosts.

Three Quantitative Trait Loci Explain More than 60% of Variation for Chill Coma Recovery Time in a Natural Population of Drosophila ananassae.

Ectothermic species such as insects are particularly vulnerable to climatic fluctuations. Nevertheless, many insects that evolved and diversified in the tropics have successfully colonized temperate regions all over the globe. To shed light on the genetic basis of cold tolerance in such species, we conducted a quantitative trait locus (QTL) mapping experiment for chill coma recovery time (CCRT) in Drosophila ananassae, a cosmopolitan species that has expanded its range from tropical to temperate regions. We created a mapping population of recombinant inbred advanced intercross lines (RIAILs) from two founder strains with diverging CCRT phenotypes. The RIAILs were phenotyped for their CCRT and, together with the founder strains, genotyped for polymorphic markers with double-digest restriction site-associated DNA (ddRAD) sequencing. Using a hierarchical mapping approach that combined standard interval mapping and a multiple-QTL model, we mapped three QTL which altogether explained 64% of the phenotypic variance. For two of the identified QTL, we found evidence of epistasis. To narrow down the list of cold tolerance candidate genes, we cross-referenced the QTL intervals with genes that we previously identified as differentially expressed in response to cold in D. ananassae, and with thermotolerance candidate genes of D. melanogaster Among the 58 differentially expressed genes that were contained within the QTL, GF15058 showed a significant interaction of the CCRT phenotype and gene expression. Further, we identified the orthologs of four D. melanogaster thermotolerance candidate genes, MtnA, klarsicht, CG5246 (D.ana/GF17132) and CG10383 (D.ana/GF14829) as candidates for cold tolerance in D. ananassae.

Direct modulation of the bone marrow mesenchymal stromal cell compartment by azacitidine enhances healthy hematopoiesis.

Mesenchymal stromal cells (MSCs) are crucial components of the bone marrow (BM) microenvironment essential for regulating self-renewal, survival, and differentiation of hematopoietic stem/progenitor cells (HSPCs) in the stem cell niche. MSCs are functionally altered in myelodysplastic syndromes (MDS) and exhibit an altered methylome compared with MSCs from healthy controls, thus contributing to disease progression. To determine whether MSCs are amenable to epigenetic therapy and if this affects their function, we examined growth, differentiation, and HSPC-supporting capacity of ex vivo-expanded MSCs from MDS patients in comparison with age-matched healthy controls after direct treatment in vitro with the hypomethylating agent azacitidine (AZA). Strikingly, we find that AZA exerts a direct effect on healthy as well as MDS-derived MSCs such that they favor support of healthy over malignant clonal HSPC expansion in coculture experiments. RNA-sequencing analyses of MSCs identified stromal networks regulated by AZA. Notably, these comprise distinct molecular pathways crucial for HSPC support, foremost extracellular matrix molecules (including collagens) and interferon pathway components. Our study demonstrates that the hypomethylating agent AZA exerts its antileukemic activity in part through a direct effect on the HSPC-supporting BM niche and provides proof of concept for the therapeutic potential of epigenetic treatment of diseased MSCs. In addition, our comprehensive data set of AZA-sensitive gene networks represents a valuable framework to guide future development of targeted epigenetic niche therapy in myeloid malignancies such as MDS and acute myeloid leukemia.

Transcriptome Analysis Reveals Candidate Genes for Cold Tolerance in Drosophila ananassae.

Coping with daily and seasonal temperature fluctuations is a key adaptive process for species to colonize temperate regions all over the globe. Over the past 18,000 years, the tropical species Drosophila ananassae expanded its home range from tropical regions in Southeast Asia to more temperate regions. Phenotypic assays of chill coma recovery time (CCRT) together with previously published population genetic data suggest that only a small number of genes underlie improved cold hardiness in the cold-adapted populations. We used high-throughput RNA sequencing to analyze differential gene expression before and after exposure to a cold shock in coldtolerant lines (those with fast chill coma recovery, CCR) and cold-sensitive lines (slow CCR) from a population originating from Bangkok, Thailand (the ancestral species range). We identified two candidate genes with a significant interaction between cold tolerance and cold shock treatment: GF14647 and GF15058. Further, our data suggest that selection for increased cold tolerance did not operate through the increased activity of heat shock proteins, but more likely through the stabilization of the actin cytoskeleton and a delayed onset of apoptosis.

Signatures of DNA Methylation across Insects Suggest Reduced DNA Methylation Levels in Holometabola.

It has been experimentally shown that DNA methylation is involved in the regulation of gene expression and the silencing of transposable element activity in eukaryotes. The variable levels of DNA methylation among different insect species indicate an evolutionarily flexible role of DNA methylation in insects, which due to a lack of comparative data is not yet well-substantiated. Here, we use computational methods to trace signatures of DNA methylation across insects by analyzing transcriptomic and genomic sequence data from all currently recognized insect orders. We conclude that: 1) a functional methylation system relying exclusively on DNA methyltransferase 1 is widespread across insects. 2) DNA methylation has potentially been lost or extremely reduced in species belonging to springtails (Collembola), flies and relatives (Diptera), and twisted-winged parasites (Strepsiptera). 3) Holometabolous insects display signs of reduced DNA methylation levels in protein-coding sequences compared with hemimetabolous insects. 4) Evolutionarily conserved insect genes associated with housekeeping functions tend to display signs of heavier DNA methylation in comparison to the genomic/transcriptomic background. With this comparative study, we provide the much needed basis for experimental and detailed comparative analyses required to gain a deeper understanding on the evolution and function of DNA methylation in insects.

Azacitidine combined with the selective FLT3 kinase inhibitor crenolanib disrupts stromal protection and inhibits expansion of residual leukemia-initiating cells in FLT3-ITD AML with concurrent epigenetic mutations.

Effectively targeting leukemia-initiating cells (LIC) in FLT3-ITD-mutated acute myeloid leukemia (AML) is crucial for cure. Tyrosine kinase inhibitors (TKI) have limited impact as single agents, failing to eradicate LIC in the bone marrow. Using primary AML samples and a patient-derived xenograft model, we investigated whether combining the FLT3-selective TKI crenolanib with the hypomethylating agent azacitidine (AZA) eliminates FLT3-ITD LIC and whether efficacy of this combination depends on co-existing mutations. Using multiparameter flow cytometry, we show FLT3-ITD occurs within the most primitive Lin-/CD33(+)/CD45dim/CD34+CD38- LIC compartment. Crenolanib alone could not target FLT3-ITD LIC in contact with niche cells while addition of AZA overcame stromal protection resulting in dramatically reduced clonogenic capacity of LIC in vitro and severely impaired engraftment in NSG mice. Strikingly, FLT3-mutated samples harboring TET2 mutations were completely resistant to crenolanib whereas neither NPM1 nor DNMT3A mutations influenced response. Conversely, primary AML LIC harboring either TET2, DNMT3A or NPM1 mutations did not show increased sensitivity to AZA. In summary, resistance of FLT3-ITD LIC to TKI depends on co-existing epigenetic mutations. However, AZA + crenolanib effectively abrogates stromal protection and inhibits survival of FLT3-ITD LIC irrespective of mutations, providing evidence for this combination as a means to suppress residual LIC.

Sex-Biased Gene Expression.

Methods of transcriptional profiling have made it possible to compare gene expression between females and males on a genome-wide scale. Such studies have revealed that sex-biased gene expression is abundant in many species, although its extent may vary greatly among tissues or developmental stages. In species with genetic sex determination, sex chromosome-specific processes, such as dosage compensation, also may influence sex-biased gene expression. Sex-biased genes, especially those with male-biased expression, often show elevated rates of both protein sequence and gene expression divergence between species, which could have a number of causes, including sexual selection, sexual antagonism, and relaxed selective constraint. Here, we review our current knowledge of sex-biased gene expression in both model and nonmodel organisms, as well as the biological and technical factors that should be considered when analyzing sex-biased expression. We also discuss current approaches to uncover the evolutionary forces that govern the evolution of sex-biased genes.

Rapid similarity search of proteins using alignments of domain arrangements.

MOTIVATION: Homology search methods are dominated by the central paradigm that sequence similarity is a proxy for common ancestry and, by extension, functional similarity. For determining sequence similarity in proteins, most widely used methods use models of sequence evolution and compare amino-acid strings in search for conserved linear stretches. Probabilistic models or sequence profiles capture the position-specific variation in an alignment of homologous sequences and can identify conserved motifs or domains. While profile-based search methods are generally more accurate than simple sequence comparison methods, they tend to be computationally more demanding. In recent years, several methods have emerged that perform protein similarity searches based on domain composition. However, few methods have considered the linear arrangements of domains when conducting similarity searches, despite strong evidence that domain order can harbour considerable functional and evolutionary signal. RESULTS: Here, we introduce an alignment scheme that uses a classical dynamic programming approach to the global alignment of domains. We illustrate that representing proteins as strings of domains (domain arrangements) and comparing these strings globally allows for a both fast and sensitive homology search. Further, we demonstrate that the presented methods complement existing methods by finding similar proteins missed by popular amino-acid-based comparison methods. AVAILABILITY: An implementation of the presented algorithms, a web-based interface as well as a command-line program for batch searching against the UniProt database can be found at http://rads.uni-muenster.de. Furthermore, we provide a JAVA API for programmatic access to domain-string­based search methods.

Quantification and functional analysis of modular protein evolution in a dense phylogenetic tree.

Modularity is a hallmark of molecular evolution. Whether considering gene regulation, the components of metabolic pathways or signaling cascades, the ability to reuse autonomous modules in different molecular contexts can expedite evolutionary innovation. Similarly, protein domains are the modules of proteins, and modular domain rearrangements can create diversity with seemingly few operations in turn allowing for swift changes to an organism's functional repertoire. Here, we assess the patterns and functional effects of modular rearrangements at high resolution. Using a well resolved and diverse group of pancrustaceans, we illustrate arrangement diversity within closely related organisms, estimate arrangement turnover frequency and establish, for the first time, branch-specific rate estimates for fusion, fission, domain addition and terminal loss. Our results show that roughly 16 new arrangements arise per million years and that between 64% and 81% of these can be explained by simple, single-step modular rearrangement events. We find evidence that the frequencies of fission and terminal deletion events increase over time, and that modular rearrangements impact all levels of the cellular signaling apparatus and thus may have strong adaptive potential. Novel arrangements that cannot be explained by simple modular rearrangements contain a significant amount of repeat domains that occur in complex patterns which we term "supra-repeats". Furthermore, these arrangements are significantly longer than those with a single-step rearrangement solution, suggesting that such arrangements may result from multi-step events. In summary, our analysis provides an integrated view and initial quantification of the patterns and functional impact of modular protein evolution in a well resolved phylogenetic tree. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.

Genomic and morphological evidence converge to resolve the enigma of Strepsiptera.

The phylogeny of insects, one of the most spectacular radiations of life on earth, has received considerable attention. However, the evolutionary roots of one intriguing group of insects, the twisted-wing parasites (Strepsiptera), remain unclear despite centuries of study and debate. Strepsiptera exhibit exceptional larval developmental features, consistent with a predicted step from direct (hemimetabolous) larval development to complete metamorphosis that could have set the stage for the spectacular radiation of metamorphic (holometabolous) insects. Here we report the sequencing of a Strepsiptera genome and show that the analysis of sequence-based genomic data (comprising more than 18 million nucleotides from nearly 4,500 genes obtained from a total of 13 insect genomes), along with genomic metacharacters, clarifies the phylogenetic origin of Strepsiptera and sheds light on the evolution of holometabolous insect development. Our results provide overwhelming support for Strepsiptera as the closest living relatives of beetles (Coleoptera). They demonstrate that the larval developmental features of Strepsiptera, reminiscent of those of hemimetabolous insects, are the result of convergence. Our analyses solve the long-standing enigma of the evolutionary roots of Strepsiptera and reveal that the holometabolous mode of insect development is more malleable than previously thought.

Inter- and intraspecific variation in Drosophila genes with sex-biased expression.

Genes with sexually dimorphic expression (sex-biased genes) often evolve rapidly and are thought to make an important contribution to reproductive isolation between species. We examined the molecular evolution of sex-biased genes in Drosophila melanogaster and D. ananassae, which represent two independent lineages within the melanogaster group. We find that strong purifying selection limits protein sequence variation within species, but that a considerable fraction of divergence between species can be attributed to positive selection. In D. melanogaster, the proportion of adaptive substitutions between species is greatest for male-biased genes and is especially high for those on the X chromosome. In contrast, male-biased genes do not show unusually high variation within or between populations. A similar pattern is seen at the level of gene expression, where sex-biased genes show high expression divergence between species, but low divergence between populations. In D. ananassae, there is no increased rate of adaptation of male-biased genes, suggesting that the type or strength of selection acting on sex-biased genes differs between lineages.

Rate of amino acid substitution is influenced by the degree and conservation of male-biased transcription over 50 myr of Drosophila evolution.

Sex-biased gene expression (i.e., the differential expression of genes between males and females) is common among sexually reproducing species. However, genes often differ in their sex-bias classification or degree of sex bias between species. There is also an unequal distribution of sex-biased genes (especially male-biased genes) between the X chromosome and the autosomes. We used whole-genome expression data and evolutionary rate estimates for two different Drosophilid lineages, melanogaster and obscura, spanning an evolutionary time scale of around 50 Myr to investigate the influence of sex-biased gene expression and chromosomal location on the rate of molecular evolution. In both lineages, the rate of protein evolution correlated positively with the male/female expression ratio. Genes with highly male-biased expression, genes expressed specifically in male reproductive tissues, and genes with conserved male-biased expression over long evolutionary time scales showed the fastest rates of evolution. An analysis of sex-biased gene evolution in both lineages revealed evidence for a "fast-X" effect in which the rate of evolution was greater for X-linked than for autosomal genes. This pattern was particularly pronounced for male-biased genes. Genes located on the obscura "neo-X" chromosome, which originated from a recent X-autosome fusion, showed rates of evolution that were intermediate between genes located on the ancestral X-chromosome and the autosomes. This suggests that the shift to X-linkage led to an increase in the rate of molecular evolution.

Dynamics and adaptive benefits of protein domain emergence and arrangements during plant genome evolution.

Plant genomes are generally very large, mostly paleopolyploid, and have numerous gene duplicates and complex genomic features such as repeats and transposable elements. Many of these features have been hypothesized to enable plants, which cannot easily escape environmental challenges, to rapidly adapt. Another mechanism, which has recently been well described as a major facilitator of rapid adaptation in bacteria, animals, and fungi but not yet for plants, is modular rearrangement of protein-coding genes. Due to the high precision of profile-based methods, rearrangements can be well captured at the protein level by characterizing the emergence, loss, and rearrangements of protein domains, their structural, functional, and evolutionary building blocks. Here, we study the dynamics of domain rearrangements and explore their adaptive benefit in 27 plant and 3 algal genomes. We use a phylogenomic approach by which we can explain the formation of 88% of all arrangements by single-step events, such as fusion, fission, and terminal loss of domains. We find many domains are lost along every lineage, but at least 500 domains are novel, that is, they are unique to green plants and emerged more or less recently. These novel domains duplicate and rearrange more readily within their genomes than ancient domains and are overproportionally involved in stress response and developmental innovations. Novel domains more often affect regulatory proteins and show a higher degree of structural disorder than ancient domains. Whereas a relatively large and well-conserved core set of single-domain proteins exists, long multi-domain arrangements tend to be species-specific. We find that duplicated genes are more often involved in rearrangements. Although fission events typically impact metabolic proteins, fusion events often create new signaling proteins essential for environmental sensing. Taken together, the high volatility of single domains and complex arrangements in plant genomes demonstrate the importance of modularity for environmental adaptability of plants.

Molecular evolution of sex-biased genes in the Drosophila ananassae subgroup.

BACKGROUND: Genes with sex-biased expression often show rapid molecular evolution between species. Previous population genetic and comparative genomic studies of Drosophila melanogaster and D. simulans revealed that male-biased genes have especially high rates of adaptive evolution. To test if this is also the case for other lineages within the melanogaster group, we investigated gene expression in D. ananassae, a species that occurs in structured populations in tropical and subtropical regions. We used custom-made microarrays and published microarray data to characterize the sex-biased expression of 129 D. ananassae genes whose D. melanogaster orthologs had been classified previously as male-biased, female-biased, or unbiased in their expression and had been studied extensively at the population-genetic level. For 43 of these genes we surveyed DNA sequence polymorphism in a natural population of D. ananassae and determined divergence to the sister species D. atripex and D. phaeopleura. RESULTS: Sex-biased expression is generally conserved between D. melanogaster and D. ananassae, with the majority of genes exhibiting the same bias in the two species. However, about one-third of the genes have either gained or lost sex-biased expression in one of the species and a small proportion of genes (approximately 4%) have changed bias from one sex to the other. The male-biased genes of D. ananassae show evidence of positive selection acting at the protein level. However, the signal of adaptive protein evolution for male-biased genes is not as strong in D. ananassae as it is in D. melanogaster and is limited to genes with conserved male-biased expression in both species. Within D. ananassae, a significant signal of adaptive evolution is also detected for female-biased and unbiased genes. CONCLUSIONS: Our findings extend previous observations of widespread adaptive protein evolution to an independent Drosophila lineage, the D. ananassae subgroup. However, the rate of adaptive evolution is not greater for male-biased genes than for female-biased or unbiased genes, which suggests that there are differences in sex-biased gene evolution between the two lineages.